N95- 22960 

ALGEBRAIC APPROACH TO THE VIBRATIONAL 
EXCITATIONS IN METHANE 

R. Lemus 1 and A. Frank 1,2 
1 Instituto de Ciencias Nucleares, UNAM 
Apdo. Postal 70-54$, Circuito Exterior, C.U. 

04510 Mexico, D.F., Mexico. 

2 fnstituto de Fisica, Laboratorio de Cuernavaca, UNAM 
Apdo. Postal 139-D, Cuernavaca, Morelos, Mexico' 



Abstract 


We present a description of the vibrational excitations of methane by means of 
an algebraic analysis of a model of coupled anharmonic oscillators. 


1 Introduction 

Consider an AB a - like tetrahedral molecule and suppose we are interested in describing its vi- 
brational degrees of freedom. This can be accomplished either in the framework of an integro- 
differential scheme or by means of an algebraic approach. The former constitutes the traditional 
method, which consists in parametrizing the Hamiltonian in terms of internal coordinates [1], 
where the potential is modeled in terms of force field constants that can be obtained from 
theoretical calculations or from fits to spectroscopic data [2]. On the other hand, the algebraic 
approach represents an alternative to the traditional methods based on the use of Lie algebras to 
represent the interactions [3]. The algebra used to describe the vibrational degrees of freedom is 
not unique. Michelot and Leroy, for example, use a unitary group U(n) as the dynamical group 
of the system with n — 1 vibrational degrees of freedom [4], while Iachello and Oss introduce 
an SU( 2) algebra for each atomic degree of freedom [5]. In this work we carry out a complete 
description of the vibrational excitations of tetrahedral molecules by assigning a 17(2) algebra 
to each interatomic potential. 


2 Algebraic Model 

The model is based on the isomorphism of the 1/(2) algebra with a one-dimensional Morse 
oscillator, whose eigenstates may be put into a one to one correspondence with a set of 17(2) D 
0(2) states, characterized by the quantum numbers |[iV],m >, as long as the value of m is 
restricted to be non-negative. In this space the Morse Hamiltonian takes the simple form Ti. — 
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•AC?20(2) j where A and N are related to the Morse potential parameters and C 2 o( 2 ) corresponds 
to the square of the 0(2) invariant operator [3]. 


B 



B 

FIG, 1. Assignment of the U l (2) algebras to tetrahedral molecules. 


For the description of a tetrahedral molecule we assign a U l { 2) algebra to each interaction 
present, as shown in Fig. 1. The first four algebras have been chosen to correspond to the A — B 
interactions, while the other six represent the B — B couplings. The molecular dynamical group 
is then given by the product U l (2) x . . . x C/ 1 0 (2), and the most general Hamiltonian, up to two 
body interactions, conserving the total number of quanta and invariant under the tetrahedral 
group 7^, can be written as 

n = n s + n B + v s ~ B . ( 1 ) 

The term H s describes the stretching degrees of freedom and has the form 


3 4 


3 4 


H s - AiY^ C 2 on 2 ) Y 0 2 O‘> (2) + Ai2 !C E &ij , 

i=l j=«+l i=l >=i+l 


z=l 


while liP is the bending contribution, given by 
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*® = *E Co* (2) + Br >% «iE E ^20^(2) + E E ^20^(2) 

j— 6 i = 5,10 j = 7,9 i=6,8 


1=5 


+ ^5,10 |^2O 510 (2) + ^20 6 8 (2) + C f 20 7 * 9 (2)| 

+ a 5 , 6 /^ y Mij+Y ^2 

(> = 6 i = 5,10 j = 7,9 i=0,8 

+ A 5i iq |At5,10 + A^6,8 + AA 7,9^ 
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The last operator, Vf _s , represents the stretching-bending interactions, which will be neglected 
as a first approximation. In these expressions C 20 * 1 ( 2 ) corresponds to the O *(2) Casimir in- 
variant, while A iij is the Majorana operator, which is related to the U l ^(2) Casimir operator 

PI. 

The simplest basis to diagonalize the Hamiltonian is the one associated to the local-mode 
chain [3] 

t/< 1 )(2) x ... x U {4 \ 2) x U { 5 ) (2) x ... x C/ (10) ( 2) D 0 ( 1 ) (2) x . . . x O ( 10 ) (2) D 0(2) 
l 11 1 i 1 1 ( 2 ) 

\[Ni] , ... [N,\ [N 2 ] , , [N 2 ] v, , , wio; V >, 

where below each group we have indicated the quantum numbers characterizing the eigenvalue 
of the corresponding invariant operator. The two boson numbers N\ and N 2 , a ^e related to the 
two sets of physical modes (stretching and bending). The quantum numbers Vi correspond to 

the number of phonons in each oscillator (1 >,■ = ^ — m*), while V = v t . 

A simple analysis of an AB 4 tetrahedral molecule [G] shows that it presents 9 vibrational 
degrees of freedom, four of them corresponding to the fundamental stretching modes (Aj 0 F 2 ) 
and the other five to the fundamental bending modes (E © F 2 ). Comparing this result with the 
local basis ( 2 ), we deduce that an unphysical bending mode is present in the algebraic formalism. 
We thus proceed to eliminate this spurious state both from the Hamiltonian and the basis. 

To accomplish this goal we first transform, for the one phonon case, the local basis to a 
normal one, which carries the irreducible representations (irreps) of the 7^ group. 'With this 
change of basis we obtain the decomposition A\ 0 F 2 for the stretches and A\ 0 E 0 F 2 for the 
bends. From this result we readily identify the Aj bending mode as the spurious state. We now 
eliminate this spurious state from the space and proceed to construct the higher phonon basis 
from the physical one-phonon set by means of the coupling coefficients C( ; ) 

v ' +v '* r y = £c(r,r 2 r; 7 i727) , ( 3 ) 

7l72 

where T and 7 label the irreps of Tj and its components, respectively. 

To eliminate the spurious contributions from the Hamiltonian we demand its null expecta- 
tion value with respect to the one-phonon spurious functions [7] 

< '•finding 1*1 ‘ ‘binding >= ° ’ (4) 


which leads to a constraint 

4(1 - N 2 )A 5 + 16(1 - 2AT 2 )B 5 ,6 + 4(1 - 2N 2 )B 5 ^ 0 = 0 
between the interaction parameters. 

The vibrational energies are obtained by diagonalizing the Hamiltonian (1) with respect to 
the normal basis ( 3 ), constructed from the projected one-phonon functions ( A \ , F 2 )— stretching 
and ( E , F 2 )- bending, taking into account the constraint (4). 
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TABLE I. Experimental [2,6,9] and calculated energies (cm 1 ) for methane. 


Normal 
V label 

r 

Expt. 

Calc. 

1 U 4 

f 2 

1310.0 

1303.7 

v 2 

E 

1533.0 

1520.4 

V\ 

A, 

2916.5 

2918.4 

^3 

F 2 

3019.4 

3027.2 

2 



2474.5 

to 

£ 

MS 


2476.4 


jRjjX 

2614.0 

2610.5 

v 2 H~ V 4 

zSI 


2827.2 


1^2 

2830.4 

2841.5 

2 v 2 



3003.7 


u 


3026.3 

"1 + "4 

f 2 

4223.0 

4222.0 


A, 


4330.9 

v 3 + ^4 

J E 


4330.9 


F , 


4330.9 



4319.0 

4330.9 

+ ^3 

(Fi 


4547.7 


1*2 

4549.0 

4547.7 

V\ + v 2 

E 


4438.8 

2 v\ 

At 


5788.0 

^1 + 

f 2 

5861.0 

5856.7 


rea 


5974.4 

2/^3 

Era 

6004.7 

6014.5 


IE 


6047.7 



4123.0 


3624.3 

3778.3 

3779.4 

3920.4 
3925.7 

3935.6 

3987.9 

4017.6 

4123.9 

4260.4 

4425.5 
4317.4 

4387.6 


Normal 

label 


I 

v\ + v 2 + 


T Expt. Calc. 


v-i + V 3 + v A 


v\ + 2v 2 

2u 2 + u 3 
2u x + v 4 
1'\ + U 3 +V A 


2v 3 + 1/4 


2 i/i + v 2 
+ ^2 + v 3 


v 2 + 2u 3 


5775.0 


5861.0 


7514.0 


5745.6 

5759.9 
5854.4 
5854.4 
5854.4 

5854.4 

5868.7 
5868.7 
5868.7 

5868.7 

5922.0 

5944.7 

6030.9 

6053.5 
6053.5 

7091.7 
7160.4 
7160.4 
7160.4 

7160.4 

7278.1 

7308.4 

7318.2 
7318.2 
7318.2 

7351.4 
7351.4 
7318.2 
7377.1 

7377.1 

7494.8 

7534.9 
7534.9 

7568.1 
7568.1 

7568.1 

8581.1 
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V 

Normal 

label 

r 

Expt. Calc. 

V 

Normal 

label 

r 

Expt. 

Calc. 



% 

4495.3 


2u\ + v 3 

f 2 

8604.0 

8603.0 


3/q < 

E 

4510.9 



A 1 


8725.5 



a 2 

4575.7 


+ 21/3 i 

f 2 

8807 

8794.1 



A, 

5392.8 



i E 


8838.5 


V\ + 2/q < 

E 

5394.7 



f 2 

8900.0 

8910.0 



1^2 

5501.7 


3^3 ^ 

F\ 


8944.8 



(Fi 

5503.6 



A ! 


8982.1 



f 2 

5503.6 



1# 

9045.0 

9034.5 



f 2 

5528.8 







tq -|- 2iq < 

f 2 

5637.7 






i 

A\ 

5637.7 








E 

5637.7 









5637.7 







3 Methane 

In this section we apply this algebraic approach to describe the vibrational levels of methane. 
According to the Hamiltonian (1) the number of parameters is eight, plus the boson numbers 
Ni and N 2 . The vibron number N\ can lie fixed from the anhannonicity of the C - H bond, 
while for the bending vibrations we have taken N 2 from the H — H interaction in H 2 O given 
in reference [ 8 ]. From these considerations, the number of free parameters is seven, taking into 
account the constraint (4). 


TABLE II. Parameters of the Hamiltonian obtained in the least square fitting 
(cm -1 ). The numbers of bosons are taken to lie N 1 =43 and N 2 =28. 


Stretching 

Bending 

rms 

A\ B 12 ^12 

— 13.2125 -0.6850 0.6328 

A 5 65,6 - 65,10 ^ 5,6 ^ 5,10 

35.4S44 2.6492 -28.0164 9.0501 5.1799 

12.16 


The Hamiltonian ( 1 ) is diagonalized in the normal basis, built by repeated couplings of 
the form (3). Since by construction this basis is symmetry adapted, the Hamiltonian matrix 
separates into blocks corresponding to the irreps of 7j. In Table I we present the least square 
fit for methane up to three quanta. Following Herzberg’s notation [ 6 ], the four fundamental 
energies for Ai, F 2 (stretching) and E, F 2 (bending) have been denoted by iq, v 2 and iq, 
respectively. The final parameters are given in Table II. The model in its simplest form (without 
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including the V B s interaction or higher order terms) seems to provide a good description of 
19 experimental energy levels with an mis deviation of 12.16 cm -1 . 


4 Conclusions 

We have preseiited a new method, which applied to an algebraic model of coupled anharmonic 
oscillators is able to describe the complete vibrational spectrum of polyatomic molecules. We 
emphasize that the method systematically incorporates group theoretical techniques which sim- 
plify the diagonalization of the Hamiltonian and provide a clear methodological procedure that 
can be applied to other molecules [10]. Although we have used the model in its simplest form, it 
can be improved in the following ways: a) Inclusion of the stretching-bending interactions V S ~ B , 
(b) Introduction in the Hamiltonian of higher order terms and (c) Addition of interactions which 
do not conserve the total number of quanta. 
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